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Abstract 


We present an analysis of the impact of heterogeneous chemistry on the partitioning of 
nitrogen species measured by the Upper Atmosphere Research Satellite (UARS) instruments. 
The UARS measurements utilized include: N 2 0, HN0 3 and C10N0 2 (Cryogen Limb Array 
Etalon Spectrometer (CLAES), version 7), temperature, methane, ozone, H 2 0, HC1, NO and 
N0 2 (HALogen Occultation Experiment (HALOE), version 18). The analysis is carried out 
for the data from January 1992 to September 1994 in the 100-1 mbar (~17-47 km) altitude 
range and over 10 degree latitude bins from 70°S to 70°N. Temporal-spatial evolution of aerosol 
surface area density (SAD) is adopted according to the Stratospheric Aerosol and Gas Ex- 
periment (SAGE) II data. A diurnal steady-state photochemical box model, constrained by 
the temperature, ozone, H 2 0, CH 4 , aerosol SAD and columns of 0 2 and 0 3 above the point 
of interest, has been used as the main tool to analyze these data. Total inorganic nitrogen 
(NO y ) is obtained by three different methods: 1) as a sum of the UARS measured NO, N0 2 , 
HNO3, and C10N0 2 ; 2) from the N 2 0-NO y correlation, 3) from the CH 4 -NO y correlation. To 
validate our current understanding of stratospheric heterogeneous chemistry for post-Pinatubo 
conditions, the model-calculated NO x /NOy ratios and the NO, N0 2 , and HN0 3 profiles are 
compared to the UARS-derived data. In general, the UARS-constrained box model captures 
the main features of nitrogen species partitioning in the post-Pinatubo years. However, the 
model underestimates the N0 2 content, particularly, in the 30-7 mbar (~ 23-32 km) range. 
Comparisons of the calculated temporal behavior of the partial columns of N0 2 and HN0 3 and 
ground based measurements at 45°S and 45°N are also presented. Our analysis indicates that 
ground-based and HALOE v.18 measurements of the N0 2 vertical columns are consistent within 
the range of their uncertainties and are systematically higher (up to 50%) than the model results 
at mid-latitudes in both hemispheres. Reasonable agreement is obtained for HN0 3 columns at 
45°S suggesting some problems with nitrogen species partitioning in the model. Outstanding 
uncertainties are discussed. 
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1. Introduction 

Prior to the launch of the Upper Atmosphere Research Satellite (UARS) in September 1991, 
nitrogen species in the stratosphere had been sampled by several different platforms, aircraft 
(e.g. recent aircraft campaigns are summarized in the following special issues: J. Geophys. 
Res., 94, 1989, nos. D9 and D14; 97, 1992, no. D8; 102, 1997, D3; Science, 261, 1128-1158, 
1993; Geophys. Res. Lett., 17, no.4, 1990; 23, no.23, 1994), balloons [e.g. Ridley et al, 1984, 
1987; Hondo et al., 1996,1997; Jucks et al, 1996; Renard et al., 1996,1997; Osterman et al., 
1997; Sen et al, 1998], ground-based [e.g. McKenzie et al., 1991; Johnston et al., 1992; Koike 
et al., 1993, 1994, 1998], limb infrared monitor of the stratosphere (LIMS) satellite [e.g. Gille 
et al, 1984; Russell et al. 1984], and the atmospheric laboratory for applications and science 
(ATLAS) series of Shuttle missions [e.g. Kaye and Miller, 1996; Rinsland et al., 1994], The 
aircraft measurements revealed many interesting features of the nitrogen species behavior in the 
stratosphere, for example, saturation of the NO*/NO y ratio as a function of aerosol surface area 
above ~5 pm 2 /cm 3 [Fahey et al, 1993] and denitrification in polar stratospheric cloud (PSC) 
regions [ Fahey et al., 1990]. However, such measurements are limited by the ER-2 ceiling altitude 
(~20 km), and their coverage is limited by aircraft paths and duration of flights. The balloon 
measurements mentioned can provide vertical profiles of several species up to 60 km [Osterman 
et al., 1997], but, again, these data have very limited coverage. The ground-based measurements 
of N0 2 and HN0 3 give the vertical profiles and column values of these species [McKenzie et 
al., 1991; Nolholt, 1994; Johnston et al., 1992; Koike et al., 1993, 1994]. Such measurements 
could provide long data sets, however, they are limited to particular stations and their analysis 
could be affected by prevailing local meteorological conditions. The LIMS sampling has given 
a near-global picture of temperature, ozone, H 2 0, N0 2 , and HN0 3 distributions during the 
October 1978-May 1979 period, which has been useful in validating global photochemical models 
[Jackman et al., 1987; Rood et al., 1993]. The ATLAS series provide the most comprehensive 
set of atmospheric species sampled, but these data have a limited temporal coverage (8-11 days) 
[Kaye and Miller , 1996]. 

Several special issues [see Geophys. Res. Lett., 19, 1992, no.2 and 20, 1993, no. 12, J. 
Geophys. Res., 98, 1993, no.D6 and 101, 1996, no.D6, J. Atmos. Sci., 52, 1995, no.17] have dis- 



cussed different aspects of the UARS measurements. Recently, Dessler et al. [1998] summarize 
principal scientific findings from the first five years of the UARS program. The UARS measure- 
ments provide a unique data set of key atmospheric species to be used for comparison of different 
satellite measurement techniques and validation of existing models of the atmosphere [Reber et 
al., 1993]. The multi-year data sets of some UARS instruments (e.g. HALOE, MLS) provide 
not only pictures of the seasonal variability of some key stratospheric species on a near-global 
basis, but also an opportunity to investigate their inter-annual changes and trends. 

Since the UARS was launched about three months after the Mt. Pinatubo eruption, there 
is an additional opportunity for atmospheric scientists to observe a near-global picture of the 
response of the atmosphere to the strongest natural perturbation in this century. Many modeling 
efforts have been undertaken to estimate different atmospheric effects of the Pintaubo aerosol 
[Boville et al., 1991; Brasseur and Granier, 1992; Pitari and Rizi, 1993; Bekki and Pyle, 1994; 
Kinnison et al., 1994; Rodriguez et al., 1994; Tie et al., 1994; Solomon et al., 1996; Rosenfield et 
al., 1997; Weisenstein et al., 1997]. Additionally, the ground-based measurements detected the 
increase of HNO3 and decrease of NO2 columns after the Pinatubo eruption [Johnston et al., 
1992; Koike et al., 1993,1994; David et al., 1994; Preston et al., 1997; Slusser et al., 1997, 1998; 
Van Roozendael et al, 1997], which could be explained by the considerably enhanced N2O5 and 
BrONC>2 hydrolysis on the Pinatubo aerosol. The HNO3 and NO2 responses to the Pinatubo 
eruption are consistent with the model predictions qualitatively. However, the magnitude of the 
observed perturbations of the HNO3 and NO2 columns at 45°S are almost two-fold stronger than 
the calculated effect by the LLNL, AER, and GSFC 2-D models [Koike et al., 1994; Rosenfield 
et al., 1997]. 

We utilize UARS data to investigate whether the observed NO x /NO y , NO, NO2, and HNO3 
distributions under the wide range of aerosol loading in the post-Pinatubo stratosphere are 
consistent with those predicted by our model calculations. We focus our study here to analysis 
of the nitrogen species only, since they play an important role for the photochemical ozone 
balance in the stratosphere. Detailed analysis of the chlorine species partitioning has been done 
by Dessler et al. [1995T996]. Our related study [Hu et al., 1998] will discuss the ozone photo- 
chemical tendencies in the stratosphere as seen by UARS. This paper shows our model analysis 
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of the UARS data acquired and incorporated in our data set. In the next section the approach 
used and model initialization are discussed. In section three we present results for the vertical 
profiles of nitrogen species, while the next section is about comparison of the column data of 
NO2 and HNO3. Finally, the last section summarizes our major findings. 

2. Approach and Model Constraints 

2.1. Approach 

A diurnal steady-state version of the AER box model without PSC chemistry is used to 
calculate the concentrations of trace species. Previous analysis of aircraft data [ Salawitch et 
al, 1994; Gao et al, 1997] shows that steady-state photochemical models can be applied almost 
everywhere in the stratosphere successfully, except in PSC processing regions or regions with 
temperature changes larger than 10-15 K/day [ Kauia et al., 1993]. A list of the included het- 
erogeneous reactions on sulfate aerosol is given in Table 1 with the updated BrON0 2 + H 2 0 
reaction probability equals to 7=l/(l/0.88 + l/exp(17. 832-0. 245wt%)) according to D.R. Han- 
son, [personal communication] (wt% is the weight percent of H 2 SC>4 in sulfate aerosol). This 
change should be noticed in the context of this paper, since the new values of 7 are larger than 
those in Hanson et al. [1996] by 0.3-0.4 and the BrON0 2 hydrolysis is not saturated even at 
high volcanic aerosol loading, thus further reducing NO x (=N0+N0 2 ) in the lower stratosphere 
[Danilin and McConnell, 1995; Lary et al., 1996]. Other reaction rates are taken from DeMore 
et al. [1997]. 

The HALOE and CLAES instruments provide 30 and approximately 1300 individual vertical 
profiles per day, respectively. To perform calculations for each profile for the period from January 
1992 to September 1994 would be an enormous task. As an alternative, we created merge files 
of the CLAES and HALOE data to represent zonal monthly or yaw-averaged conditions (see 
section 2.2). These files are used to provide input for the model runs. The measured species from 
these files are compared against the model results. To check the validity of this approach, we 
performed daily runs for February 1993. We picked this month because atmospheric circulation 
is perturbed at the end of the winter in the northern hemisphere, thus providing potentially a 
worst case scenario for our study. We calculate the monthly averaged NO x /NO y ratios using 





our merge files as input. The value of <^§*1 daii y > averaged over the month is compared with 


a single calculation for the monthly averaged <7^§ 1 >|mont/i/y at the same level and latitude bin. 


Figure 1 depicts the relative difference of this comparison for 
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measurements, respectively. In general, the difference between daily and monthly runs is in 


the ±10% range almost everywhere, except for some points in the lower stratosphere. We ana- 
lyzed the most striking points where the difference was above the ±10% range. Our analysis of 
these points (not shown here) reveals that the individual daily NO x /NO y ratios could be very 
scattered. For example, among 15 individual HALOE measurements on February 22 and 23 at 
15°N and 68 mb, 4 measurements give NO x /NO y ~ 0.25, 7 measurements between 0.17 and 
0 .2, 4 measurements show very low NO x /NO y <0.05 (thus giving < jf^\daiiy >=0.16), while 
< Wc? > |montwy = 0.2. The other points outside the 10% difference range in Figure 1 are also 
characterized by a large scattering of individual HALOE measurements. Our conclusion from 
this comparison is that: 1) our approach using monthly averaged merge files is reasonable for 
the purpose of this study and 2) analysis of the individual UARS profiles, if necessary, should 


be done very cautiously, especially in the lower stratosphere, where the accuracy of individual 
measurements could be low (15% for HALOE NO and NO2, 6-15% for CLAES HNO3, 23% for 
CLAES CIONO2 [Morris et al., 1997; Table 1]) and less data may be available. 


2.2. UARS Merged Data 

We created a merged UARS data set containing the zonally-averaged stratospheric trace 
species on daily, monthly and yaw-averaged (i.e. during the 36-day yaw cycle) bases. Daily 
CLAES v.7 measurements released on CD-ROMs and the recently released v.18 of the HALOE 
data are used to create these merge files. The files contain measured concentrations of the 
following species: HALOE CH 4 , NO, N0 2 , H 2 0, 0 3 , HC1, and temperature; CLAES N 2 0, 
HNO3, CIONO2, and N2O5; at the following 13 pressure levels: 100, 68.2, 46.4, 31.6, 21.5, 
14.7, 10, 6.8, 4.6, 3.2, 2.2, 1.5, and 1 mbar. The daily files simply summarize available HALOE 
measurements andCLAES data within ±5° latitude and ±10° longitude distance from a HALOE 
point during a particular day at these levels. 
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There are 10 yaw-averaged files per year, since the duration of a yaw cycle is approximately 
36 days. We bin the data by latitude with a 10° step from 70°S to 70°N and mark the UARS 
results by a mean latitude value in each bin (i.e. 65°S corresponds to the 70°S-60°S band 
etc.). Available UARS data beyond the 70°S-70°N range (namely, 80°S-70°S and 70°N-80°N) 
are ignored because of the following two reasons: first, they are rather sparse and cannot provide 
enough coincident measurements by the CLAES and HALOE instruments in a 10° bin due to 
their sampling geometry; second, our diurnal steady-state model approach may be less valid in 
these regions, especially after polar vortex breakups. The HALOE measurements for each yaw 
cycle and latitude bin are averaged, thus forming a skeleton of the merge files for the sunrise 
and sunset data, separately. Only those CLAES daily measurements that are within 5° latitude 
distance from a HALOE sampling on the same day are selected for the yaw-averaged files. They 
are included in the merge files after simply zonally averaging. This approach should be robust 
for the long-lived species, like HN0 3 and N 2 0. Possible errors in C10N0 2 data (its diurnal 
variability could contribute only 3% at most in NO y [Moms et al., 1997]) have a negligible 
effect on the conclusions of this study. Performing a similar analysis, Morris et al. [1997] 
concluded that a more sophisticated method (like trajectory mapping technique [Morris et al., 
1995]) is not necessary for such kinds of studies, since the zonal averaged data provide almost the 
same results. Each yaw-averaged merge file also contains variances of the measured parameters 
during a yaw cycle. 

The same approach is used to create monthly zonally averaged merge files. Additionally, we 
created files of zonally monthly averaged CLAES measurements of N 2 0, HN0 3 , and C10N0 2 , 
disregarding their proximity to the HALOE points. These files of nitrous oxide zonally averaged 
data <N 2 O c > are used to calculate NO y , Cl y , and Br y to initialize our box model runs for the 
N 2 0 approach. Our analysis is focused below on the zonally averaged merge files created for 
each yaw cycle from January 1992 to September 1994 and related to the sunset and sunrise 
HALOE measurements. 

2.3. Model Constraints 

We initialize our calculations with H 2 0, CH 4 , and 0 3 data from the yaw-averaged merge 
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files described above. Ozone is kept constant in the model and equal to its HALOE value. Our 
merge files contain also SAGE II aerosol SAD, ozone and oxygen columns above each level, and 
averaged julian day of the HALOE measurements to calculate photolysis rates. To calculate 
diurnal behavior of the solar zenith angle, solar declination is adopted for the mean day of 
the HALOE measurements for the yaw cycle and latitude bin considered. The ozone column 
above the point of interest is obtained by integrating the available HALOE ozone profile for 
that particular latitude and month. The O 2 column is defined by the pressure level of interest 
p as O^ olumn (cm -2 ) = 4.444xl0 21 xp (in mbar). 

Since we are running our box model until steady-state conditions, initial partitionings within 
NO y , Cl y and Br y families are not important. However, UARS does not measure all nitrogen, 
chlorine and bromine species required to constrain total inorganic nitrogen (NOy), inorganic 
chlorine (Cly), and inorganic bromine (Br y ) for the model runs, so methods of model initializa- 
tion are needed. Three alternative methods to constrain these families are outlined in Table 2 
and discussed in the following subsections. 

2.3a. NOy 

The first way is to define total nitrogen NOy (= NO -I- NO2 + NO3 -I- HNO3 -I- 2xN 2 0s + 
HONO -I- HN0 4 + CIONO2 + BrON0 2 ) as NO^ = < NO u > + < NO^ > + < HNO3 > + 
< ClONOj >, where the superscript U denotes UARS measurements of NO, N0 2 (by HALOE), 
HNO3, and C10N0 2 (by CLAES), and the symbol <> denotes zonally averaged value of the 
bracketed species during a yaw cycle. We did not include CLAES N2O5 in the equation for 
NOy (since these measurements are recommended for use only at 3.2, 2.2, and 1.5 mbar [Kumer 
et al., 1996b]), thus underestimating NO y in this approach. The other omitted species in the 
equation for NO^ (NO3, HONO, HNO4, BrONOz) may contribute only a few percent or less to 
NOy, and thus are unimportant for the purposes of this study. 

The second approach is to use zonally averaged values of CLAES <N 2 O c > to derive NOy 
using equation (1) in the northern hemisphere by Loewenstein et al ., (1993] and equation (2) in 
the southern hemisphere by Keim et al. [1997]: 

NO y = 20.7- 0.0644 x N 2 0, (1) 

NO y = 21.82 - 0.0699 x N 2 0, (2) 
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in the above equation NO y and N 2 0 are in ppbv. These empirical N 2 0-NO y correlations are 
valid for N 2 O>50 ppbv (below 10 mbar) [ Fahey et al ., 1990; Prather and Remsberg, 1993]. 
Above this altitude, destruction of NO y via photolysis of NO followed by the N + NO > N 2 
+ O reaction cause the correlation curve to turn over. Above these altitudes, we use NO y 
or utilize the correlations derived from ATMOS (see Figure 2 and discussion below). These 
two approaches can be used to initialize our box model only until May 1993, since the CLAES 
cryogen supply was exhausted after that time. 

The third method, which uses the HALOE measurements of CH 4 , allow us to extend our 
model analysis after May 1993. This approach uses the empirical correlation between NO y , 
N 2 0, and CH 4 based on the analysis of ATLAS/ATMOS data [ Michelsen et al., 1998]. Their 
analysis provides two sets of the equations: 

in the tropics (i.e. between 20°S and 20°N) for 0.4 ppmv < CH 4 < 1.28 ppmv: 

NO„ = -1.3519xl0“ 8 + 0.05794(CH 4 ) - 18544(CH 4 ) 2 - 7.3339xl0 9 (CH 4 ) 3 , 
for CH 4 > 1.28 ppmv: 

NO y = -7.9107xl0~ 7 + 1.4051(CH 4 ) + 8.3839xl0 5 (CH 4 ) 2 - 1.6787xl0 n (CH 4 ) 3 , 
and outside tropics for CH 4 <0.9 ppmv: 

NO y = -1.2321 x 10 -9 + 0.02276(CH 4 )- 31866(CH 4 ) 2 - 4.178xl0 10 (CH 4 ) 3 ; 
for CH 4 >0.9 ppmv: 

NO y = 9.8251 xl(T 9 + 0.0418(CH 4 )- 55088(CH 4 ) 2 + 1.608xl0 10 (CH 4 ) 3 , 

here NO y and CH 4 are in mixing ratio. We must stress the approximate nature of these cor- 
relations. They do not account for possible latitudinal and seasonal changes, which could be 
particularly important at higher altitudes [ Keim et al., 1997]. There are several situations when 
the HALOE NO* measured at higher altitudes (where NO x is the major NO y component) is 
larger than the NO y derived from the CH 4 correlations. In these cases, the initial NO y is as- 
signed to be equal to NO"' ALOB + HNO" LAES + C10N0" LAE5 or to N0" ALOE (if CLAES data 
are not available). 

Figure 2 compares the NO y fields in January 1993 calculated from the three methods out- 
lined above. This month is chosen for illustration because it contains the largest latitudinal 
coverage by the sunset UARS measurements. The top panel shows NO" (= HALOE sunset 
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monthly zonally averaged values of <NO> and <N0 2 > together with CLAES <HN0 3 > and 
<C10N0 2 >). The middle panel depicts the difference between NO y , inferred from the zonal 
yaw-averaged N 2 0 field by CLAES <N 2 O c >, and NOy . The bottom panel shows the difference 
between NO y , obtained according to the above CH 4 -NO y correlation for the HALOE <CH 4 >, 
and NO^ . In general, all three approaches have similar features: relatively low NO y values in the 
lower stratosphere, especially in the tropics; gradual increase with altitude up to the maximum 
values above 20 mbar and decrease above 3 mbar. However, more detailed analysis also reveals 
principal differences in these NO y fields. For example, the CH 4 approach (bottom panel) gives 
the highest level of NO y with its maximum above 18 ppbv (i.e. 1-2 ppbv higher than NO y ) in 
the tropics at 10-3 mbar (~34-40km). Compared to NOy , the CH 4 approach overestimates the 
NO y content up to 3 ppbv in the upper stratosphere in the tropics, near the tropical barrier 
and winter mid-latitudes up to 3 ppbv and underestimates up to 3 ppbv above 2 mbar in the 
summer mid-latitudes. Switching from tropical to extratropical CH 4 -NO y correlations, we see 
some problems near 20°S and 20°N (lower panel). The N 2 0 approach (middle panel) shows, in 
general, a good agreement with the UARS approach. However, somewhat larger NO y mixing 
ratios (>1 ppbv) are obtained below 40 mbar at southern mid-latitudes due to low values of 
CLAES N 2 0 there. Also, the N 2 0 method provides smaller NO y values than those in the UARS 
approach by >1 ppbv (middle panel) between 10-7 mbar in the tropics. Above 7 mbar, the NO y 
contents are the same for the N 2 0 and UARS approaches. Obviously, the UARS approach has 
the smallest altitudinal-latitudinal coverage below 30 mbar (since it requires almost coincident 
measurements of all four species: NO, N0 2 , HNO3, and C10N0 2 , while the other methods 
require only N 2 0 CL/4E5 and CH ^ ALOE ), thus potentially leaving more gaps in the model analy- 
sis. In general, the three approaches give reasonably close results at mid-latitudes. The largest 
discrepances are observed in the tropics and subtropics and in the upper stratosphere, maybe 
due to a combination of interference from Pinatubo aerosols and poorly defined tracer-tracer 
correlations. 

Our analysis of similar plots for other months and years shows that the above mentioned 
differences and similarities of the NO y field are persistent. For example, figure 3 show the 
difference between NO y distribution during the first yaw cycles (mostly occurred in January) 
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in 1992 and 1993 defined by the UARS (top panel), CLAES N 2 0 (middle panel), and HALOE 
CH 4 (bottom panel) methods. This figure depicts the change in atmospheric dynamics from 
January 1992 to January 1993 (since NO y could serve as a tracer) and compares three methods 
to define NOy. There are the following similarities: higher values of NOy in January 1992 below 
10-20 mbar (30-40 mbar) in the southern (northern) mid-latitudes; higher NOy content above 2 
mbar in the tropics in January 1992. Among the differences, one can mention stronger change of 
NO y below 40 mbar at northern and southern mid-latitudes and different values of NOy change 
between 30°S and 40°N at the pressure range of 30-5 mbar in the CH 4 scheme. 

Since the UARS does not measure NOy itself, Figures 2-3 illustrate only the quality of 
the above correlations and, strictly speaking, cannot conclude which approach is right. In 
general, NOy derived directly from UARS is preferable in the lowermost stratosphere, provided 
that' the error bars for individual components are not too large. The UARS initialization of 
NOy is preferable in the upper stratosphere, since the N 2 0-NOy and CH 4 -NO y correlations 
are less reliable there. Initialization of NO y from the output of a global model for a particular 
location, season, and altitude could introduce additional uncertainties due to model formulations 
of dynamics and chemistry. For example, Nevison et al. [1997] reported that the Garcia-Solomon 
model overestimates NOy in the upper stratosphere. 

Our sensitivity study shows that the model NO x /NOy ratio is insensitive to NO y , changing 
roughly a couple percent per 10-20% of NO y change. This helps us to ignore possible uncer- 
tainties of the NOy initialization, analyzing the NO x /NO y behavior. 

2.3b. Cly and Br y 

Initial inorganic chlorine Cl y (=HC1 + C10N0 2 + HOC1+ Cl + CIO + BrCl + 2(C1 2 
+ C1 2 0 2 )) is defined from the correlation of total (organic and inorganic) chlorine with N 2 0 
[ Woodbndge et al., 1995; their equation (11)] and dependence of the Cl y fraction of total chlorine 
as a function of N 2 0 [their Figure 7]. The initial values of inorganic bromine Br y (= Br + BrO 
+ HOBr 4- BrON0 2 + BrCl + HBr + 2Br 2 ) are inferred from the empirical correlation of Br y 
with N 2 0 [ R.J. Salauntch, personal communication, 1997]: 


11 



for 120 ppbv< N 2 0 < 312 ppbv: 


Br y = 13.009 + 0.1284(N 2 O) - 5.452 x 1(T 4 (N 2 0) 2 , (3) ' 

for N 2 0 < 120 ppbv: 

Br y = 21.0984 - 0.5317(N 2 0/120) (4) 

where Br y is in pptv and N 2 0 is in ppbv (Br y =l pptv is assumed at N 2 0> 312 ppbv). These 
values of Br y could differ by up to 2 pptv from the recent analysis by Wamsley et al. [1998] in 
the stratosphere. However, this potential difference is not important for the analysis of nitrogen 
species presented here. To calculate the Cl y and Br y contents, the N 2 0 values could be obtained 
from the yaw-averaged CLAES <N 2 O c > data before May 1993 or HALOE data of <CH 4 > 
after May 1993 using the CH 4 -N 2 0 correlation [Michelsen et al., 1998]: 
in the tropics (i.e. 20°S-20°N): 

N 2 0* = -1.2851 x 10“ 7 + 0.7742(CH 4 )- 1.5510x 10 6 (CH 4 ) 2 + 1.3097xlO 12 (CH 4 ) 3 - 3.3888x 10 17 (CH 4 ) 4 
outside the tropics: 

N 2 0* =-2.2309 xl0~ 8 + 0.1934(CH 4 )- 5.6278 x 10 5 (CH 4 ) 2 + 7.3019xl0 n (CH 4 ) 3 - 2.34xl0 17 (CH 4 ) 4 , 
here N 2 0* and CH 4 are in mixing ratio. The contents of inorganic chlorine and bromine species 
could be different for the two approaches, when both CLAES N 2 0 and HALOE CH 4 are avail- 
able. However, the possible differences do not significantly affect nitrogen species partitioning 
investigated in this paper. 

3. Analysis of Vertical Profiles of Nitrogen Species 
3.1. NOx/NOj, 

We performed model runs for each yaw period from January 1992 to September 1994 for 
sunrise and sunset measurements, except for June 1992, when HALOE did not operate because 
of technical problems aboard the UARS. Latitudinal coverage of our merge files and, conse- 
quently, our model calculations differs from month to month, depending on the geometry of 
the HALOE and CLAES samplings. In this study we focus on analysis of the nitrogen species 
behavior primarily at 45°S and 45°N, since additional ground-based measurements of N0 2 and 
HN0 3 are available at these latitudes. Such a comparison of model calculations with measure- 


12 



ments obtained by two very different techniques should be a very useful test to analyze the 
representation of stratospheric chemistry in the models. 

Figure 4 shows temporal evolution of the sunset NO x /NO y ratio and aerosol SAD at four 
levels: 68, 46, 31, and 22 mbar levels at 45°S (left column) and 45°N (right column). We show 
these levels because the lower stratosphere has the strongest perturbation due to Pinatubo 
aerosol [e.g. Yue et al, 1994; Lambert et al . , 1997; Thomason et al . , 1997], Results at the 
lowermost level of 100 mbar are not shown because the UARS measurements are very noisy 
due to retrieval problems within a dense Pinatubo aerosol cloud. Since the model NO x /NO y 
ratio is insensitive to the NO y values, our model results are shown by red crosses in Figure 4 for 
the CH 4 approach only. However, since the UARS does not measure NO y itself (while NO x is 
measured by the HALOE instrument), there are differences between UARS NO x /NO y results 
because of the different ways to initialize NO y according to the UARS (cyan symbols), N 2 0 
(blue symbols), or CH 4 (green symbols) approaches. The principal differences between these 
initializations are already discussed in the previous section. The monthly variable deviations 
shown are derived from statistical analysis of the available UARS data at each level during a 
particular yaw period. 

A seasonal modulation of the NO x /NO y ratio is clearly seen in both UARS data and model 
results. This ratio is maximum during summer (January at 45°S, June- July at 45°N) and 
minimum during winter (July at 45°S, January at 45°N) because of the increased photolysis of 
HN0 3 in summer and the increased heterogeneous processing via N 2 0 5 hydrolysis in winter. 
Since there are many gaps in the UARS data set at this latitude, it is hard to see an interannual 
increase of the NO x /NO y ratio due to the shown reduction of SAD after the Pinatubo eruption. 
However, the NO x recovery can be seen more clearly in the model results. Missing UARS points 
at 68 mbar during the first half of 1992 in Figure 4 are due to deterioration of the HALOE NO x 
measurements by the Pinatubo aerosol at this level [ Gordley et al., 1996]. 

Comparison of the model and UARS results shows a good agreement at the 68 and 46 mbar 
levels throughout the whole period analyzed. However, at the upper levels (particularly, at 22 
mbar), a clear underestimation of the NO x /NO y ratio by the box model in summer in both 
hemispheres is evident. At some points (January 1993 at 45°S and July 1992 at 45°N, 22 mbar) 
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this disagreement can be more than 30%. A similar problem was also noticed by Morris et al 
[1997; Plate 1), who applied the Goddard 2-D model with the UARS-constrained photochemistry 
and concluded that “the model NO x /NO y ratio is marginally but systematically lower than the 
UARS ratio, especially at 800K (~31 km) and 1000K (~35 km)”. Recent analysis of the satellite 
[e.g., Lary et al., 1997] and balloon measurements [Sen et al., 1998; Kondo et al., 1998] also 
show that current photochemical models with heterogeneous chemistry underestimate the NO x 
amount in the stratosphere. 

To investigate the altitudinal-latitudinal dependence of the NO x /NO y ratio, we present re- 
sults in Figure 5 for January 1993 between 68 and 4.6 mbar. The UARS NO x /NO y ratio is shown 
by the cyan, blue and green symbols for the UARS, N 2 0, and CH 4 approaches, respectively 
(see Table 1). The model values of NO x /NO y , calculated for the CH 4 approach, are shown by - 
red and black lines for calculations with and without heterogeneous reactions on sulfate aerosol, 
respectively. Increase of the NO x /NO y ratio towards the summer hemisphere at all levels is due 
to faster photolysis of HNO3 there. We performed model runs with gas phase chemistry only to 
investigate a sensitivity of the model NO x /NO y ratio to the heterogeneous reactions on sulfate 
aerosol. Recently, Donahue et al. [1997] argue that current recommendations [DeMore et al., 
1997] overestimate the HNO3 formation reaction rates by 10-30%. Additionally, we run our 
model with heterogeneous reactions on sulfate aerosol and with the reduced by 25% reaction 
rate of OH -i- N0 2 + M — > HNO3 + M (shown by the grey lines) to reflect the possibility 
outlined by Donahue et al. [1997], assuming that their data are valid for the stratospheric 
conditions. 

The agreement between model calculations with heterogeneous chemistry and UARS mea- 
surements is reasonable below 31 mbar, especially in the winter (i.e. northern) hemisphere. In 
the summer (i.e. southern) hemisphere our model tends to underpredict the NO x content at and 
above 46 mb. The gas phase chemistry results (black lines) clearly overestimate the NO x /NO y 
ratio at these levels. However, at three levels, 21.5, 14.7 and 10.0 mbar, model results with 
heterogeneous chemistry underestimate the NO x /NO y ratio by 0.1-0. 2 at all latitudes. Even 
the NO x /NO y results for the gas phase chemistry runs fall below the measurements at 10 and 
14.7 mbar everywhere and at 21.5 mbar at southern mid-latitudes. Since the aerosol SAD drops 
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with altitude, the difference between model results with and without heterogeneous chemistry 
is diminished at higher levels (black and red curves are coincident at 6.8 and 4.6 mb). At the 
two upper layers, 4.6 and 6.8 mbar, model results are consistent with the UARS measurements. 
The 25% reduction of the OH + N0 2 4- M — ■» HN0 3 + M helps to increase the NO./NOj, 
ratio. However, this increase is still too small to explain the UARS measurements. 

In our opinion, the most alarming signal for modelers from Figures 4 and 5 is a striking 
underprediction the NO x /NOj, ratio in the lower and middle stratosphere, since predictive ca- 
pabilities of global models could be affected, particularly, for tasks where accurate knowledge of 
NO and N0 2 at these altitudes are important or perturbations of interest include NO x emissions 

(like aircraft). 

3.2. NO, N0 2 , HN0 3 

To investigate model underestimation of the NO x /NO y ratio in the stratosphere more care- 
fully, we analyze separately the sunset vertical profiles of NO, N0 2 , and HN0 3 in January 1993 
at 45°S (summer) and 35°N (winter) according to the UARS initialization (see Figure 6). The 
UARS values are shown by symbols, the box model profiles for the UARS approach by solid 
and dashed lines for calculations with and without heterogeneous reactions on sulfate aerosol, 
respectively. Below 30 mbar (~ 24 km) and above 6 mbar (~ 34 km), the model reproduces 
the observed profiles of these nitrogen species reasonably well. Results for N 2 0 5 and C10N0 2 
are not shown here, since 1) the maximum N 2 Os mixing ratio at sunset is small (a few tenth of 
ppbv) due to its diurnal variability and 2) the C10N0 2 profile has a relatively small maximum 
(~ 1 ppbv at 20 mbar) which is properly reproduced by our model. 

However, the box model considerably underestimates N0 2 and overestimates the HN0 3 
amount in the 30-7 mbar range by several ppbv. The reason for this model underestimation is 
not clear yet and further analysis is required. Additional missing heterogeneous mechanism of 
renoxification on sulfate aerosol, like HN0 3 4- H 2 CO — > 2HONO + C0 2 + H 2 0 [Fairbrother 
et al., 1997; Iraci and Tolbert, 1997] could barely help, since the sulfate aerosol SAD is small 
at these levels (e.g., <0.1 fF /cm 3 at 10 mbar in January- 1993 at 45°S). Possible renoxification 
reactions on soot [Tabor et al., 1994; Rogaski et al., 1997; Lary et al., 1997; Bekki, 1997] also 
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could not help, since there is not a noticeable amount of soot at these altitudes, according to 
current knowledge. Since results for the gas phase run (black lines) fall below the UARS data at 
10 and 14.7 mbar, the heterogeneous chemistry could hardly be responsible for this disagreement. 
Sen et al [1998] also reported that the balloon measurements of NO2 exceed the values calculated 
by the steady-state model by 20-30% at 15-30 km and this difference is larger than the estimated 
total uncertainty of the measurements. This disagreement could be reduced (but not solved!) by 
10% changes of some photochemical constants (e.g. NO+O3 — > N0 2 + 0 2 and N0 2 absorption 
cross section) within a given range of their uncertainties by [DeMore et al,, 1997]. Recent 
POLARIS campaign results show that photochemical models underpredict the NO x amount in 
the lower stratosphere during summer under almost background aerosol loading [R.J. Salamtch, 
personal communication , 1998]. It is worth noticing that this model shortcoming in the altitude 
region between 30 and 7 mbar is responsible for consequent disagreement between models, 
on the one hand, and UARS and ground-base measurements, on the other hand, of the N0 2 
columns (see next section). Perhaps something is missing in our current understanding of the 
gas phase or photolytic processes with nitrogen species at these levels. A careful look at the 
reaction rates crucial for the NO x and HN0 3 partitioning is required. Additionally, Zip f and 
Prasad [1998] speculated that an oxygen molecule above 25 km could get excited by a photon in 
the Schumann-Runge bands and react in this excited short-lived state with N 2 producing two 
molecules of NO. 

Since the HNO3 content in the stratosphere is expected to increase after the volcanic erup- 
tions due to enhanced N 2 0 5 hydrolysis [ Rodriguez et al., 1991, 1994; Brasseur and Granier, 
1992; Bekki and Pyle , 1994; Kinnison et al, 1994; Risland et al., 1994], it is interesting to see 
how the CLAES HN0 3 distribution did react to the Pinatubo eruption to check for consistency 
between UARS and the model picture of this event. Figure 7 shows the change of the HN0 3 
distribution between January 1992 and January 1993 as seen by the CLAES instrument (top 
panel) and as predicted by our model with heterogeneous reactions on sulfate aerosol for the 
UARS initialization (bottom panel). The corresponding change in the NOj, distribution for the 
different model initializations is shown in Figure 2. Figure 7 depicts the difference of the nitric 
acid distribution obtained by simple subtraction of the CLAES zonally monthly averaged HN0 3 
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profiles in January 1993 from the similar values in January 1992. A similar procedure is applied 
to the model calculated HNO3 fields during the same period. These changes reflect a complex 
interplay between dynamics and chemistry of the stratosphere from January 1992 to January 
1993. Below we make an attempt to separate effects of dynamics and chemistry in the HNO3 
field. 

The pattern of the nitric acid change shows a strong interhemispheric difference. The area 
of the largest increase of HNO3 (more than 2 ppbv) as seen by the CLAES instrument is 
located in the southern hemisphere between 35 and 10 mbar. Our UARS-initialized box model 
(bottom panel of Figure 7) captures its location and magnitude very well, albeit the area of 
the calculated HN0 3 change above 2 ppbv is smaller. The heterogeneous chemistry on sulfate 
aerosol (primarily due to N 2 0 5 + H 2 0 — > 2HNO3) is responsible for the HN0 3 change in this 
area, since the change in NO y (which could serve as an indicator of a change of stratospheric 
dynamics) is smaller (top panel of Figure 2) than the change in HNO3. Also, CLAES and 
model calculations show positive change in the HNO3 field almost everywhere in the southern 
hemisphere stratosphere, while negative change of the NO y field is seen there above 10 mbar. 
In the northern mid-latitudes, positive change of HNO3 distribution is coincident with negative 
change of the NO y field. Since the opposite signs of the HNO3 and NO y changes at the same 
location imply strong repartitioning of the nitrogen species from January 1992 to January 1993 
due to heterogeneous chemistry, we can conclude that the whole stratosphere has been strongly 
chemically perturbed during this period of time. Above 4 mbar the changes of HNO3 (both 
CLAES and model) and NO y are small. 

Analysis of the change in the NO y distribution above 20 mbar in the southern mid-latitudes 
shows a very weak interannual change of dynamics there from 1992 to 1993 (see Figure 2). 
On the other hand, the situation in the northern mid-latitudes is quite different, perhaps, due 
to the asymmetric QBO-like effects [ Gray and Ruth, 1992; Kumer et al., 1996a]. Changes in 
NO y larger than -4 ppbv (inferred from CLAES N 2 0) are present in a quite noticeable part of 
the stratosphere above 50 mbar in northern tropics and mid-latitudes. However, the HALOE 
instrument has seen somewhat smaller perturbation in the CH 4 field in this region as shown in 
the bottom panel of Figure 2 via the CH 4 -NO y correlation. 
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Model NOj, results for the N 2 0 and CH 4 approaches (Figure 2, middle and bottom panels) 
show some similarities with Figure 7. However, they also demonstrate that dynamical pertur- 
bation of the long-lived tracers such as N 2 0 and CH 4 (translated into NOj, change in the middle 
and bottom panel of Figure 2) was so strong in some locations (especially near 20 mbar at 20°N) 
that its fingerprint was observed in the HNO3 field as a negative change (Figure 7) overshadow 
effects of heterogeneous chemistry there. 

4 . Discussion of the HNO 3 and N0 2 Column Data 

Ground-based concurrent measurements of HN0 3 and N0 2 columns at Lauder (New Zealand, 
45°S, 170°E) provide an additional opportunity to validate our model calculations and UARS 
measurements. Figure 8 compares monthly averaged ground based measurements (blue symbols 
with errorbars) of HN0 3 (top panel) and N0 2 columns at sunset (bottom panel) with UARS 
data and our box model calculations with the three different initializations. The CLAES daily 
measurements of HN0 3 at 44°S (version 8) and HALOE sunset monthly values of N0 2 between 
50° and 40°S are shown by green. The CLAES instrument does not account for nitric acid below 
100 mbar. To compare CLAES data with ground-based observations, we increased the partial 
HN0 3 column values by 0.9 xlO 15 cm -2 as in Kumer et al [1996a], assuming ~70 pptv HN0 3 
below 100 mbar. Following McKenzie et al. [1991], we assume a constant tropospheric N0 2 
profile of 20 pptv in summer and 60 pptv in winter at 45° [ Logan et al., 1981] with sinusoidal 
inter-seasonal change, thus adding a value between ~2.8xl0 14 and ~7.7xl0 14 cm -2 , respec- 
tively, to the HALOE N0 2 column values. The same values of the HNO3 and N0 2 columns 
below 100 mbar have been added to our box model calculations. 

The ground-based N0 2 column measurements were made using a visible spectrometer based 
on the twilight zenith-sky technique. The N0 2 slant column amount was measured at sunrise and 
sunset at solar zenith angle of 90 degrees. The vertical column value was calculated assuming an 
air mass factor of 16.7. The absolute accuracy estimated only from the temperature dependence 
of the N0 2 absorption cross-section was -1-18 ± 5 % (overestimate), because the 20°C values 
were used for the analyses. Uncertainty in the air mass factor was estimated to be ± 12 %. 
Other sources of the systematic error are described in Koike et al. [1998]. Vertical column HNO3 
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amount was measured using Fourier transform infrared spectrometer. The absolute accuracy 
estimated from uncertainties in the HNO3 absorption line parameter and vertical HNO3 profile 
was ± 13 % [Jones et al, 1994], 

The ground-based and CLAES measurements of the HN0 3 agree well from January 1992 to 
May 1993, while the box-model calculations are consistent with the measurements after August 

1992 for all initializations. The CH 4 - and N 2 0-initialized box model overestimates the vertical 
HNO3 amount in the first half of 1992 due to overestimate of the NO y content below 30 mbar 
inferred from the N 2 0 and CH 4 fields and shown in Figure 1. Our UARS initialization can 
provide the column value only during the first (^January) and third («March) yaw cycles in 

1993 (shown by black symbol), primarily, because of the missing HALOE N0 2 measurements 
at the lowermost levels (100 and 68 mbar) affected by the Pinatubo cloud. Both ground- 
based and CLAES measurements capture a noticeable increase of the HN0 3 column in summer 
1992 at the top of its seasonal cycle due to enhanced NO^ conversion on Pinatubo aerosol. 
Unfortunately, our box model could not be initialized during May-July 1992 due to missing 
HALOE measurements at this latitude. 

Analysis of the N0 2 column data (bottom panel) demonstrates a recovery of N0 2 after the 
Pinatubo eruption and good agreement between ground-based and HALOE sunset measure- 
ments. Our box model calculations, as could be anticipated from our section 3 and previous 
comparisons [Koike et al., 1994; Rosenfeld et al., 1997], considerably underestimates the N0 2 
column up to 50% with the ground-based measurements during the whole period of time in- 
vestigated. As we show in Figure 6, the reason of this discrepancy is mainly due to the model 
underestimate of the N0 2 content between 30 and 6 mbar. Assuming 50% N0 2 underestima- 
tion in the model in the 20-35 km altitude range (which can contributes ~75% of the total N0 2 
column [McKenzie et al, 1991; their Figure 8] in summer, we increase our model column values 
by ~38%, making them much closer to the ground-based and HALOE measurements. The new 
version 18 of HALOE data is in better agreement with the ground-based measurements than 
the version 17 results [Gordley et al., 1996]. 

Results for the sunrise N0 2 column show a very similar pattern with lower absolute values 
by a factor of ~1.6 and are not shown here. Satisfactory agreement of the HN0 3 column 
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values and disagreement of N0 2 column values between model calculations and measurements 
indicates problems in the model’s ability to reproduce proper partitioning between NO x and 
HN0 3 rather than an incorrect initialization of NO y . A similar problem was obtained in the 
AER 2-D model calculations for the post-Pinatubo conditions (not shown here). Also, our 
box model showed very similar results as the GSFC 3-D model and both models underpredict 
the NO content in the stratosphere during analysis of the balloon measurements in September 
1994 [Kondo et al, 1998], As we mentioned above, several other models with state-of-the-art 
heterogeneous chemistry and updated photochemical kinetics underestimated the NO x content 
in the stratosphere [e.g. Lary et al., 1996; Morris et al., 1997; Sen et al, 1998]. 

Figure 9 shows a comparison of the ground-based (Moshiri station, Japan, 44°N), HALOE 
and model values of the N0 2 column at 45°N which depicts similar problems outlined in our 
discussion for 45°S. The HALOE v.18 N0 2 column agrees well with the ground-based measure- 
ments. Our model calculations are smaller than the ground-based measurements by about 50% 
through the whole year. A similar pattern is obtained for sunrise conditions. The reasons for 
these discrepances are outlined above. 

5. Conclusions 

We have presented analyses of the UARS measurements (HALOE and CLAES) using the 
AER diurnal steady-state box model from January 1992 to September 1994. It was shown that 
use of zonal yaw-averaged UARS measurements, which avoid high uncertainties of the UARS 
individual profiles and considerably reduce computer time for analysis, captures the principal 
features of the nitrogen species partitioning in the global post-Pinatubo stratosphere. 

The AER diurnal steady-state photochemical stratospheric model, with proper constraints 
from the UARS measurements, can be rather successfully applied for analysis of this data set. 
We discussed three possibilities listed in Table 1 to initialize the box model from the available 
set of UARS measurements. Each of them has particular advantages and disadvantages. Our 
analysis shows that the UARS and N 2 0 approaches can provide reasonable agreement in the 
NOx/NO,, ratio between model calculations and UARS data (especially in the upper strato- 
sphere). However, they are limited by the CLAES operation time (until May 1993). The CH 4 
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approach has a larger temporal coverage, but is has several outlined weaknesses, like overesti- 
mation of NO y above 10 mbar in the tropics. 

In general, our results show that the model and the approach used provide good agreement 
with the UARS data below 30 mbar, thus supporting our current parameterization of the het- 
erogeneous reactions on sulfate aerosol. Our analysis, concentrated at 45°S and 45°N, revealed a 
seasonally modulated gradual increase of the NO x /NOy ratio after the Pinatubo eruption in the 
lower and middle stratosphere, consistent with the UARS measurements. Our state-of-the-art 
module of heterogeneous chemistry on sulfate aerosol reproduces the vertical profiles of NO, 
N0 2 , and HN0 3 along with the temporal evolution of the latitudinal-altitudinal distribution of 
these species. 

Among clear discrepances between model results and UARS measurements, we emphasize 
a clear underestimation of the N0 2 content at the 30-7 mbar levels (particularly evident in the 
summer hemisphere) in our box model. Since the aerosol loading is very low at these altitudes, 
it signals possible problems in gas-phase chemistry or photolytic processes. This problem has 
been noted in previous studies [e.g. Considine et al., 1992; Moms et al., 1997; Sen et al, 
1998], however, there is no clear understanding yet of why it occurs. Analysis of the N0 2 
column measurements show that the ground-based (at Lauder, New Zealand, 45°S and Moshiri, 
Japan, 44°N) and HALOE v.18 data are consistent within the range of their uncertainties, 
while the model calculations underestimate the N0 2 amount by up to 50% during the period 
and latitudes considered. Further analysis is urgently needed to understand the reason for the 
model underestimation of N0 2 at these altitudes (~24-34 km) and its consequencies for global 
2-D (3-D) model predictions. 
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Table 1. Heterogeneous reactions on sulfate aerosol included in the model. The reaction proba- 
bilities are taken from DeMore et al. [1997], except the BrON0 2 + H 2 0 reaction (see text). 

N 2 0 5 + H 2 0 — > HN0 3 + HN0 3 

BrON0 2 + H 2 0 — > HOBr + HNO 3 
CIONO 2 + H 2 0 — > HOC1 + HN0 3 

C10N0 2 + HC1 — » Cl 2 + HN0 3 
HOC1 + HC 1 — > Cl 2 + H 2 0 
HOBr + HC1 — > BrCl -I- H 2 0 


Table 2 . Initialization methods for NOy, Cly, and Br y employed in this study, using the 
UARS data. 


SPECIES 


APPROACH 



UARS 

N 2 0 CLi4£5 

CH? ALOE 

NOy below 7 mbar 

NO" 

fi(N 2 O c ) 

f 4 (CH 4 ) 

NO y above 7 mbar 

NO" 

NO" 

f 4 (CH 4 ) 

Cly 

f 2 (N 2 O c ) 

f 2 (N 2 O c ) 

f 2 (N 2 0*) 

Br y 

f 3 (N 2 O c ) 

f 3 (N 2 O c ) 

f 3 (N 2 0*) 


NO^ =<NO t/ >+<NO^ >+<HNO^ >+<C10N0^ >, N 2 O c is the monthly CLAES N 2 0 data, 
N 2 0* is the N 2 0 value from the CH 4 -N 2 0 correlation (see text); function fi(N 2 0) is defined by 
Keim et al [1997] and by Loewenstein et al. [1993] in the southern and northern hemisphere, 
respectively; functions f 2 (N 2 0), f 3 (N 2 0), and f 4 (CH 4 ) are defined according to Woodbridge et 
al. [1995], R.J. Salawitch, and Michelsen et al. [1998], respectively (see text). 

Figure Captions 

Figure 1 . Difference (in %) between model NO x /NO y calculated for monthly zonally averaged 
UARS data and for daily UARS data averaged over month and longitude. The CH 4 approach 
to initialize model runs has been used. A gap between 25°S and 5°N is due to the geometry 
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of HALOE measurements in February 1993 (sunrise and sunset measurements in the southern 
and northern hemisphere, respectively). 

Figure 2. Distribution of NO y (in ppbv, top panel) in January 1993 using the UARS initialization 
(see Table 1). The middle and bottom panels show the difference between UOy ARS and NO y 
as inferred from the correlations with CLAES zonal monthly averaged <N20 C > and HALOE 
<CH 4 >, respectively (see text for details). 

Figure 3. Change of NO y between January 1992 and January 1993 inferred from the different 
methods of NO y initialization according to: the UARS approach (top panel); the N 2 0 approach 
(middle panel), and the CH 4 approach (bottom panel). 

Figure 4. Temporal evolution of the sunset NO x /NO y ratio from January 1992 to September 1994 
at 45°S (left column) and 45°N (right column) at the 68, 46, 32, and 22 mbar levels. Green, 
cyan, and blue symbols with errorbars correspond to the UARS data with NO y according to the 
CH 4 , UARS, and N 2 0 initializations, respectively; red crosses to box model calculations. The 
N 2 0 and UARS approaches can provide NO y data only until May 1993 because of the CLAES 
lifetime. The black lines show behavior of aerosol surface area density (right vertical axis, in 
/im 2 /cm 3 ). 

Figure 5. Latitudinal dependence of the sunset NO z /NO y ratio in January 1993 at the levels 
shown between 68.2 and 4.6 mbar. Model results with and without heterogeneous reactions on 
sulfate aerosol are shown by red and black lines, respectively. Model results with heterogeneous 
chemistry and with reduced by 25% reaction rate of OH + N0 2 + M — * HNO3 + M are 
shown by grey lines. The UARS values are shown by symbols for the three methods of NO y 
initialization: cyan for NO y (UARS), green for NO y (CH 4 ), and blue for NO y (N 2 0). 

Figure 6. Vertical profiles of NO (blue), N0 2 (cyan), and HN0 3 (red) in January 1993 at 45°S 
(top panel) and 35°N (bottom panel) at sunset. UARS data are shown by symbols with er- 
rorbars, model results are shown by solid (with heterogeneous reactions) and dashed (without 
heterogeneous reactions) lines according to the UARS initialization. 

Figure 7. Change of the zonally averaged HNO3 distribution (in ppbv) from January 1992 to 
January 1993 according to the CLAES measurements (top panel) and model calculations for 
the UARS initialization (bottom panel). 
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Figure 8. HNO3 (top panel) and sunset NO2 (bottom panel) columns at 45°S between Jan- 
uary 1992 and January 1995. Blue symbols show ground-based measurements in Lauder (New 
Zealand); red, pink and black symbols- box model results for the CH 4 , N 2 0 and UARS ini- 
tializations, respectively. Green curves in the top panel depict daily HNO3 column values at 
44°S according to the CLAES instrument (version 8). Green symbols in the bottom panel show 
monthly- zonally averaged HALOE (v.18) values of the sunset N0 2 column. The box model 
and UARS data are increased by partial HNO3 and N0 2 columns below 100 mbar (see text). 

Figure 9. Comparison of the N0 2 column (in 10 15 cm~ 2 ) at sunset according to the ground-based 
measurements at Moshiri (44°N, blue), our box model calculations (red - CH 4 approach, pink 
crosses- N 2 0 approach, and black- UARS approach), and the HALOE monthly-zonally averaged 
data (green) from January 1992 to January 1995. The sunrise measurements (not shown) depict 
similar pattern with ~ 1.6 times smaller values. 
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